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ABSTRACT 

The morphology, kinematics and entrainment mechanism of the HH 46/47 molecular outflow were 
studied using new ALMA Cycle observations. Results show that the blue and red lobes are strik- 
ingly different. We argue that these differences are partly due to contrasting ambient densities that 
result in different wind components having a distinct eflFect on the entrained gas in each lobe. A 
29-point mosaic, covering the two lobes at an angular resolution of about 3'^ detected outflow emis- 
sion at much higher velocities than previous observations, resulting in significantly higher estimates 
of the outflow momentum and kinetic energy than previous studies of this source, using the CO (1-0) 
line. The morphology and the kinematics of the gas in the blue lobe are consistent with models of 
outflow entrainment by a wide-angle wind, and a simple model may describe the observed structures 
in the position- velocity diagram and the integrated intensity map. The red lobe exhibits a more 
complex structure, and there is evidence that this lobe is entrained by a wide-angle wind and a colli- 
mated episodic wind. Three major clumps along the outflow axis show velocity distribution consistent 
with prompt entrainment by different bow shocks formed by periodic mass ejection episodes which 
take place every few hundred years. Position-velocity cuts perpendicular to the outflow cavity show 
gradients where the velocity increases towards the outflow axis, inconsistent with outflow rotation. 
Additionally, we find evidence for the existence of a small outflow driven by a binary companion. 
Subject headings: ISM: jets and outflows — stars: formation — ISM: Herbig-Haro objects — ISM: 
individual (HH 46, HH 47) 



1. INTRODUCTION 

As stars form inside molecular clouds, they eject mass 
in energetic bipolar outflows. The resulting bipolar wind 
from a young stellar object (YSO) may reveal itself 
through Herbig-Haro (HH) objects observed in the opti- 
cal, H2 emission knots in the infrared (IR), and molecular 
(CO) outflows observed at millimeter (mm) wavelengths. 
HH objects delineate highly collimated jets and their (in- 
ternal or leading) bow shocks. The H2 IR emission also 
arises from recently shocked gas and in many cases it 
traces the bow-shock wings that extend towards the driv- 



ing source. CO outflows map the ambient gas that has 
been swept-up well after it has been entrained by the 
protostellar wind and has cooled. Hence, these differ- 
ent manifestations provide complementary views: while 
HH objects and H2 emission provide a "snapshot" of the 
current shock interaction, the CO outflow trace th e pro- 
tostar's mass loss history (e.g.. Richer et al.||20QQ ). 

Protostellar winds inject energy and mornentum into 
the surroundings, thereby perturbing the star-formation 
environment. Outflows may be responsible for the cle ar- 
ing of material from the core (Arce & Sargent 2006), a 



Arce, et al. 



process that could result in the termina tion of the infall 
phase (e.g., [Velusamy fc Langer||1998L affect the star 
formation etiicien cy in the cloud (e.g., Matzner fc Mc- 
iKee 2000; Nakam ura"XTil[2QQ7| [Machida & Hosokawa 
p013), and determine the niass of stars (Myers 2008). In 
addition, outflows can affect the kinematics, density and 
chemistry of a substantial volume of their parent clouds, 
and thus can be important to the turbulent dynamics and 
energetics of their host cores. How exactly protostellar 
winds entrain and disperse the surrounding gas and feed 
their parent clouds turbulence remains a mystery. High 
angular (and velocity) resolution, mult i- wavelength, ob- 
servations of outflows are needed to probe all different 
outflow manifestations in order to compare them with 
hydrodynamic numerical simulations and shock models 
to obtain a complete picture of protostellar winds, the 
entrainment process and their impact on their surround- 
ings. 

Here we present Atacama Large Millimeter/sub- 
millimeter Array (ALMA) observations of the HH 46/47 
molecular outflow. The high resolution data were ob- 
tained with the aim of studying the entrainment process, 
the underlying protostellar wind and the properties of the 
molecular outflow. The observed outflow is driven by HH 
47 IRS (a.k.a., HH 46 IRS 1, IRAS 08242 - 5050), a low- 
mass Class I YSO with a total luminosity of about 12 L© 
in the Bok globule ESO 2 16-6 A, located on the outskirts 
of the Gum Nebula at a distance of 450 pc ISchwartz 



surrounding the source show the molec ular outflow gas 



1977! "Stanke et an[T999! 'Reipurth et~aL]|2QQQ; Noriega- 



5^ 



r^respo et al.|2004p. The source lies very close to the edge 



of the Bok globule, which explains why the blue (north- 
ern) lobe of the HH 46/47 bipolar flow is clearly seen 
at optical wavelengths outside the globule, while within 
^ 2' of the source the red lobe (which mostly lies inside 
the globule) is best seen in IR images. Hubble Space 
Telescope (HST) observations indicate that HH 47 IRS 
is actually a binary system, where the two components 
of the system are observed to be separa ted by only 0.^^26 
or about 120 AU (iReipur th et~ar]|2000[ ). 

Many of the properties of the HH 46/47 flow have 
been determined through extensive optical and IR ob- 
servations. The combination of optical spectral data and 
proper motion studies of the HH knots (using ground- 
based telescopes) allowed an estimate of the flow's incli- 
nation to the plane of the sky (^ 30°) and an average 
jet velocity of 300 km s~^ ( Reipurth 1989; Reipurth 
Heathcote||199T| 



^^^^... ..^^^^JEisloeffel fc Mundt..l994, Micono et al 

1998p . More recently, proper motion studies (using HST) 



of the blue lobe optical HH knots, in combination with 
the r esults from sp ectrosco pic Fabry-Perot observations 
(^Morseet al."1994), allowed 'Hartigan et al^ ([2QQ5) to es- 
timate an average orientation angle, with respect to the 
plane of the sky, of about 37° . Infrared images show how 
shocked H2 emission in the southern (redshifted) lobe 
traces the walls of a 2' long and 36^^ wide cavity pro- 



duced by the flow (Eisloffel et al. 1994 Noriega-Crespo 
|et al. 2004). 

The HH 46/47 molecular outflow has only been fully 
mapped at low resolution (with beams larger than 24'') 
by several a uthors, usi ng the t wo lowest rotational transi- 
tions of CO ([C hernin fc Massonll991[ [Olberg et al.|1992[ 
npen et al. 2009|). AFEA observations of higher 
6) of a smaller area (80'' x 80") 



van Kempen 

(JO transitions (up to 
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reaches temperatures as high as 100 K (van Kempen et 
ar]|200 9'). Moreover, recent Herschel PACS observations 
of a region within 30" from the source reveal the exis- 
tence of hot CO and H2O (thought to be partly produced 
by non-dissociative C-shocks) and OH and [01] emission 
from dissociative J-shocks where the protostellar wind 
interacts with the surround ing de nse core ( |Van Kempen] 
et al.||2010[|Wampfler et al.||2010D . 

Even though the extensive multi-wavelength studies of 
this source have helped understand the physics of proto- 
stellar jets, the lack of high-resolution molecular outflow 
observations (until now) has hampered our ability to ob- 
tain a complete picture of the outflow phenomenon in 
HH 46/47, the entrainment process and the protostellar 
wind's impact on the cloud. Here we present the first in- 
ter ferometric, high-resolution (beam ~ 3") observations 
of the HH 46/47 molecular (CO) outflow and compare 
the morphology and kinematics to existing (simple) mod- 
els of outflow entrainment. In a future paper we will use 
hydrodynamical simulations to model the outflow and 
aim to provide more stringent constraints on the launch- 
ing and entrainment mechanisms in this flow. 

2. OBSERVATIONS 

The observations were carried out using ALMA from 
28 December 2011 to 25 January 2012, during the Early 
Science Cycle phase. The Band 3 data were obtained 
over eleven scheduling blocks, with 16 to 18 antennas in 
the (Cycle 0) compact configuration, consisting of pro- 
jected baselines in the range of 12 to 277 m. The corre- 
lator was configured to observe three windows with the 
highest spectral resolution {^v = 30.5 kHz), with a band- 
width of 58.6 MHz and centered at the 115.27 GHz, the 
frequency of the CO (1-0) transition, and at 100.88 and 
100.07 GHz. We only detected line emission in the 115.27 
GHz window, which has a velocity resolution of 0.08 km 
s~^. No line was detected in the other two windows, 
and were therefore used to map the continuum emission. 
We mapped our sources using a rectangular 29-point mo- 
saic, with contiguous pointings separated by 24.5", ori- 
ented at a position angle of about 58°, designed to cover 
the length of the HH 46/47 outflow from HH 47A, in the 
northeast, to HH 46C, in the southwest, and the width of 
the red outflow lob e seen in IR images (see, e.g., Noriega- 
Crespo et al.||20 04'). The mosaic had three rows, the two 
outer rows consisted of 10 pointings and the central row 
had 9, resulting in a map of about 1.2' x 4.2', centered 
on 08^25^41^5, -51°00'47" (J2000). 

The nearby quasar J0845-5458 was used for phase and 
gain calibration. Flux calibration was carried out us- 
ing observations of Mars, and J0538-440 and 3c279 were 
used as bandpass calibrators. The visibility data were 
edited, calibrated and imaged (using the CLEAN algo- 
rithm) in CASA, using 0.5" cells and 1024 cells in each 
spacial dimension. For the spectral data we deflned a dif- 
ferent clean region for each channel, encircling the area 
with the brightest emission. The resulting synthesized 
beam for the ^^00(1-0) data cube is 3.2" x 2.4", with 
a PA= 16.5°, and the (la) rms is 25 mJy beam~^ in a 
(spectroscopically smoothed) map with a velocity reso- 
lution of 0.5 km s~^. The continuum was obtained by 
averaging the emission over all channels in the two win- 
dows close to 100 GHz with no detected line emission, 
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Fig. 1. — Integrated intensity map of the HH 46/47 CO(l-O) outflow. Blue (red) contours represent the blueshifted (redshifted) lobe. The 
blue lobe is integrated over —6.8 < Vout < —1.3 km s"-*^. The lowest contour and subsequent contour steps are 0.2 and 0.25 Jy beam"-*^ km 
s"-*^, respectively. The red lobe is integrated over 1.3 < Vout < 6.7 km s"-*^. The lowest contour and subsequent contour steps are 0.25 and 
0.3 Jy beam"-*^ km s"-*^, respectively. Green contours show the 3 mm continuum emission surrounding HH 47 IRS. The lowest contour and 
subsequent contour steps are 1.3 and 1 mJy beam"-*^, respectively. The thick grey line delineates the edge of the globule as traced by our 
CO data close to the cloud velocity (see Figure |3|. The synthesized beam of the CO map is shown on the lower left corner of the figure. 
Grey diagonal parallel dashed lines show the edges of our map. The position of redshifted outfiow clumps Rl, R2 and R3 is shown, as well 
as that of the blueshifted clur np BO and clump A. The dark dashed line passing through BO, HH 47 IRS and A indicates the position of 
the p — V cut shown in Figure [iT] 



resulting in a bandwidth of 117.2 MHz. The synthesized 
beam and rms for the continuum map are 3.1'^ x 2.2^' 
(PA= 15°), and 0.33 mJy beam"! 



3. RESULTS 

3.1. Continuum Source 

We detected continuum thermal dust emission at 100 
GHz near the position of the protostar, with a peak in- 
tensity of 8.1 mJy beam~^. A Gaussian fit puts the con- 
tinuum peak at 8^25^43.^8,-51°00'36'' (J2000.0), which 
we adopt as the position of HH 47 IRS, and it is con- 
sistent with that given by Reipurth et al. ( 2000| ). The 
continuum emission is barely resolved, and it is slightly 
extended toward the southwest (se e Figure T]) . Th e total 
flux density is 14 mJy. Following ISchnee et al.l ([2Q10|) 

Kempei 



and assuming a dust temperature of 30 K (|van Kemp en 



let al. 2009), we obtain a total (gas and dust) mass aS' 
sociated with the continuum emission of 0.9 Mq. The 
angular resolution of our continuum is not enough to re- 
solve different peaks associated with the two different 
binary components of the system. The observed con- 
tinuum emission most probably traces the envelope that 
surrounds the binary system, and observations with an 
order of magnitude higher angular resolution are needed 
to investigate whether each component has its own sep- 
arate circumstellar envelope. 



3.2. Line Data and the Molecular Outflow 

Integrated intensity maps of the CO (1-0) blueshifted 
and redshifted emission from the HH 46/47 molecular 
outflow are displayed in Figure [l] The most striking as- 
pect of this outflow is the clear difference in size and mor- 
phology between the two different lobes; the red (south- 
western) lobe extends approximately 2' from the source, 
while the blue (northeastern) lobe only extends up to 
about 30^' from the powering YSO. This asymmetry had 
been noted in earlier, lower angular resolution, observa- 
tions of the CO outflow and it is generally assumed that 
it is due to the fact that the H H 47 IRS is close to the 
edge of the parent globule (i.e, [C hernin fc Mass on||1991 
piberg et aL]|1992[ jvan Kempen et al..2009} . The proto- 
stellar wind's blueshifted (northeastern) lobe breaks-out 
of the cloud, where there is little molecular gas for it to 
entrain, while the redshifted lobe dives into the globule 
and is able to entrain much more gas along its path. 

In addition to the drastic difference in length, each 
lobe exhibits dissimilar morphologies most likely caused 
by differing dominant entrainment processes. The blue 
lobe shows mainly a parabolic morphology, with a posi- 
tion angle of about 60° . The morphology and the velocity 
structure of the blue lobe suggests this molecular outflow 
is mostly being formed by the entrainment of cloud ma- 
terial by a wide angle protostellar wind (see § |4.1|). We 



also detect in the blue lobe a small protuberance south 
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Fig . 2. — Comparison between CO(l-O) spectrum from |01berg| 
|et al.||T992| ) and our ALMA observations. The spectrum shown m 
green is the same one shown in Figure 4 of Olberg et al. (1992), 
which comes from the brightest position in the redshifted lobe of 
their the map. The black histogram shows the spectrum at the 
same position from a map using our ALMA data that has been 
smoothed to the same beam size as that of the Olberg et al. data 
(44''). The spectrum from the ALMA data is also smoothed to a 
velocity resolution of 0.5 km s"-*^ for easier visual comparison with 
the Olberg et al. spectrum. 

east of the continuum source (labeled BO in Figure [l]) , 
which we argueis due to the outflow from a binary com- 
panion (see § 4.4). The red lobe, on the other hand, has a 
V-shape near the source, with an opening angle of about 
65°. Starting at about 30^' from the source, the inte- 
grated intensity contours mostly trace the southern and 
northern outflow cavity walls, which extend up to 80'^ 
and 2', respectively, from the source. We suspect that 
these structures trace the limb-brightened walls of an 
approximately cylindrical shell produce by the outflow- 
cloud interaction. In addition, we detect three distinct 
clumps of gas along the outflow axis which we label Rl, 
R2 and R3 (see Figure [l]). In § 4.2 we argue that these 
result from approximately periodic mass ejections driven 
by the central protostar. 

In Figure [2] we plot the spectrum at the position of 
brighte st integrated intens ity in the CO(l-O) outflow 
map of Olberg et aL (|1992|). We comp are the spectrum 
shown in i^lgure 4 of Olberg et al. ( 1992[ ) with a spectrum 
obtained from our CO map, smoothed to the same reso- 
lution as the beam size of the Olberg et al. observations 
(44^^). From this comparison we see that the ALMA ob- 
servations recover most (if not all) of the emission at LSR 
velocities greater than 6.5 km s~^ and less than 3.5 km 
s~^, and it is clear that we do not fully recover the emis- 
sion at velocities close to that of the cloud {vlsr ^ 5 km 
s~^, van Kempen et al. 2009). We are therefore confident 
that our ALMA data can be used to reliably study the 
emission at outflow velocities greater than about 1.5 km 
s-i. 

3.2.1. Channel Maps 

Figure |3] shows the emission over the outflow velocity 
range between approximately —10 and 10 km s~^, and 
reveals how the CO emission structure of the region de- 
pends on velocity. We define outflow velocity (vout) as 
the LSR velocity of the emission minus the cloud LSR 
velocity. Panels in the upper row (blueshifted velocities) 
and lower row (redshifted velocities) present the emission 
averaged over a 1 km s~^ range centered on the velocity 



indicated in the upper right corner of each panel. Panels 
in the center row show emission averaged over 0.5 km 
s~^ wide range for velocities close to that of the cloud. 
We note that low intensity outflow emission is detected 
at higher outflow velocities than those shown in Figure [3J 
This weak high- velocity emission is only seen when inte- 
grating the emission over a larger velocity range (> 5 km 



4.2 



s "j, see § |4.1[ and 

In Figure |5| the first" two panels in the upper row 
i^^out = —5.8 and —4.8 km s~^) show an approximately 
symmetric parabolic (low-intensity) structure extending 
northeast of the source. At lower blue outflow velocities 
the emission is stronger and the northern wall of the lobe 
extends further from the source than the southern wall. 
At Vout = — 2.8 km s~^ we detect two faint arcs of emis- 
sion, about 40'^ northeast and 2' southwest of the proto- 
star, which we ascribe to the limb-brightened edge of the 
host globule. These arcs are brighter closer to the cloud 
velocity (see panels at Vout = —0.1 and 0.4 km s~^). At 
^out = — 1.8 km s~^ the blueshifted outflow emission is 
very bright and it is concentrated in an inclined V-shape 
structure northeast of the source. It is only at very low 
outflow velocities that we detect the clump BO, about 
20^' southeast of HH 47 IRS, which very likely traces the 
molecular outflow driven by the binary compon ent that 
is not driving the main HH 46/47 outflow (see § |4.4[ ). In 
this panel we also detect faint emission southwest of the 
outflow source, coincident with some of the regions where 
we detect redshifted outflow emission. It is possible that 
the walls of the southwest (mainly redshifted) outflow 
lobe are expanding into the surrounding cloud (similar 
to the RNO 91 outflow, Lee & Ho 2005), and the fila- 
mentary blueshifted emission comes from the front walls 
of the cavity. In addition, at Vout = —1.8 km s~^ we 
detect a small relatively bright clump close the position 
of the HH 47C shock (and the redshifted clump R3) at 
about 2' southwest from the source, which most likely 
arises from gas entrained by HH 47C. 

As shown in Figure [2J close to the cloud velocity our 
interferometer observations do not fully recover the emis- 
sion from the cloud, as the observations are not sensitive 
to large scale (> 20^') structures. This is the reason why 
the channel maps from Vout = —1.1 to 0.9 km s~^ appear 
so clumpy and the emission at these velocities cannot be 
used to obtain reliable estimates of the outflow or cloud 
mass nor to map the overall distribution of the cloud 
molecular gas. There are, however, a few small-scale fea- 
tures that are discernible among the clumpy mess. Even 
at Vout = —1-1 km s~^ we are able to identify the V- 
like structure just north-east of HH 47 IRS due to the 
blue outflow lobe, and a clump of gas close the position 
of HH 47C. At the velocity of the cloud (i.e., the Vout 
= —0.1 and 0.4 km s~^ panels) we detect emission from 
the edge of the cloud, which is coincident with the diffuse 
8 jam. emission observed in Spitzer Space Telescope im- 
ages of_the_£egi^n__di^^^_bord^r^ edge of the host glob- 
ule dNoriega-Crespo et aL]|2004[ fV^lusamy et al.,,2007^ . 
This structure seems to trace the outer parts of the cloud 
that are being heated by the UV radiation from the stars 
that produce the Gum Nebula. The Vout = 0-4 km s~^ 
panel also exhibits a bright structure that extends north- 
east of the so urce and coincides with the HH 46 optical 
nebula (see § 4.1). This vey low velocity emission could 
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3.2.2. Outflow Mass 

To obtain a reliable estimate of the outflow mass we 
would need to know the opacity of the ^^ CO (1-0) line, 
as i n outflows the opacity is expected to vary with veloc- 
ity dBally et al.||1999[ [Yu et aL|[T999l |Arce fc Goodman 
2001b). One way to estimate the ^^CO line opacity is us- 



8''25"'48' 45' 42' 39' 36' 33' 
RA (J2000) 

Fig. 3. — Channel maps of the CO(l-O) emission. Top (bottom) row shows blueshifted (redshifted) emission, where each panel shows 
the emission average over 1 km s"-*^. The lowest contour and subsequent contour steps are 0.075 and 0.15 Jy beam"-*^, respectively. The 
middle row shows velocities near the cloud velocity, and each panel shows the emission average over 0.5 km s"-*^. The lowest contour and 
subsequent contour steps are 0.15 and 0.3 Jy beam"-*^, respectively. The central outflow velocity is given in the upper right corner of each 
panel in units of km s"-*^. The synthesized beam is shown on the lower right corner of each panel. The red cross shows the position of HH 
47 IRS. 

be due to the gas that is entrained by the expanding 
back wall of the northeast (mainly blueshifted) lobe. At 
Vout = 0.9 km s~^ the overall structure of the CO emis- 
sion is similar to the integrated intensity image of the 
redshifted lobe in Figure [T] yet there are some notable 
differences. These include a bright narrow structure that 
extends northeast of the source (coincident with the blue 
lobe) , and the fact that the wide angle morphology of the 
southwest lobe is more prominent here t han i n the inte- 
grated intensity map (which we argue in 
for the existence of a wide-angle wind). 

Similar to the blue lobe, the redishfted outflow lobe 
is brighter at lower outflow velocities. The Vout= 1-7 
km s~^ panel shows the strongest emission, where the 
highest intensity is seen close to the source and on the 
southern and northern walls of the outflow lobe. In this 
panel the CO emission shows a wide-angle morphology 
close to the sources, and the emission associated with the 
northern outflow wall is somewhat extended (as opposed 
to a simple narrow structure as seen at higher outflow 
velocities). The three redshifted clumps (Rl, R2 and 
R3) are seen most clearly at Vout =2.7 and 3.7 km s~^, 
but they are also detected out to Vout ^ 5 km s~^. The 
channel maps also show that from Vout =1.7 km s~^ to 
"^out = 5.7 km s~^ the redshifted lobe becomes slightly 
narrower as the velocity increases. 



4.2 is evidence 



ing the line ratio of ^^CO to ^^CO. Unfortunately, there 
are no existing observations of the ^^ CO (1-0) emission 
with an angular resolution (and sensitivity) similar to our 
^^CO map that would allow us to properly correct for the 
opacity. We thus estimate the mass assuming the line is 
optically thin, and warn that our mass estimates may 
be underestimated, from a fact or of few up to as much 
magnitude (e.g., [Cab rit & Bertout] 
ian'l2001bl lOflner et al. 2011). We : 



as an order of 



1990 
follow 



Arce & Goodman 
Bourke et al. (19971 
from the CO (1 

[^^C0]/[H2] = 10 , . „ . , . ^ 

that if we consid er the same velocity ranges used by 01- 



to estimate the outflow properties 
"ata, a ssuming an abundan ce ratio of 



as m 



Olberg et al. (1992). We note 



berg et al. ( |1992 ) for the two different outflow lobes (i.e. 
—5.0 to —1.5 km s~^ for the blue lobe and 1.5 to 7.5 km 
s~^ for the red lobe), we obtain a total outflow mass of 
3.6 X 10~^ Mq, similar to the value of 0.03 M© obtained 
by Olberg et al. This is consistent with Figure [2J which 
shows that our interferometric observations recover most 
(if not all) of the emission at outflow velocities greater 
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TABLE 1 

Outflows Properties 




Lobe 


Tex = 


Mass^ [10-2 M©] 
15 K 50 K 100 K 


Momentum^' ^ [10-^ M© km s"^] 
15 K 50 K 100 K 


Energy^' ^ [10^^ ^^.g] 
15 K 50 K 100 K 


Blue 
Red 




12 3 

3 8 14 


3/6 8/16 15/30 
12 / 24 27 / 54 51 / 102 


3/12 8/32 14/56 
9/36 21/84 39 / 156 



^ Estimates are obtained assuming emission is optically thin. Values should be treated as lower limits. 

Values before the slash are not corrected for the outflow inclination and those after the slash are corrected assuming 
an inclination of the outflow axis, with respect to the plane of the sky, of 30° . 



than 1.5 km s~^. 

An estimate of the CO outflow excitation temperature 
is needed to calculate the outflow mass. However, there 
is a very wide range of temperature values giv en in the 
literature for the HH 46/47 molecular outflow. Chernin 



fc Masson] ( 1991 ) use the ratio of brightness temperature 
between tne (JU(3 — 2) and CO (2 — 1) at outflow veloci- 
ties to estimate an outflow excitation temperature {Tex) 
of 8.5 ±1 K (but use a value of 10 K in their outflow mass 
calculations). While Olberg et al. (11992) estimate an ex- 
citation temperature, at outflow velocities, of 15 K using 
the in tensity ratio of CO (1 - 0) to CO (2-1). On the other 
hand, van Kemp en et ah] ( |2009[ ) estimate much larger 
values tor '1'^^ (between 50 and 150 K, but most likely 
about 100 K) using the CO(3-2)/CO(6-5) ratio. With- 
out ALMA observations of higher CO transition lines 
we cannot make our own estimate of Tg^;, we therefore 
calculate the outflow mass, and other outflow proper- 
ties shown in Table [T] using different values of excitation 
temperature. It is clear that uncertainties in T^x result 
in highly inaccurate estimates of the outflow properties. 
Our observations are much more sensitive than any 
previous study of the HH 46/47 molecular outflow, which 
allows us to detect outflow emission over a much wider 
velocity range. In this work we measure the outflow mass 
(yMout), momentum [Pout = ^ Mout{vout) Vout], and en- 
ergy [Eout = 0.5 S Mout{vout) vl^t] over the velocity 
range where we detect outflow emission with intensity 
greater than 3cr, that is —30 < Vout < —1.6 km s~^ 
for the blue lobe and 1.4 < Vout < 40 km s~^ for the 
red lobe. In both lobes the very high velocity outflow 
(I'^on^l > 20 km s~^) is constrained to within about 15'^ 
of the source. Detection of low-intensity CO (1-0) outflow 
emission at very high velocities resulted in our estimates 
of the outflow kinetic energ y an d momentum to be signif- 
icantly higher than those of Olberg et alj(il992|. This im- 
plies that other similar molecular outflows may be much 
more energetic than what previous (low-sensitivity) ob- 
servations indicated and outflows may have the potential 
to have more impact on their surrounding cloud than pre- 
viously thought. 

4. DISCUSSION 

4.1. Blue lobe: evidence for wide-angle wind 

As mentioned above, we detect blueshifted outflow 
emission out to larger velocities than what is presented 
in Figure [3J The intensity of the high- velocity emission is 
very low and it is necessary to sum over velocity ranges 
wider than 1 km s~^ in order to detect significant emis- 
sion. In Figure [4] we show four velocity maps, integrated 
over different velocity ranges, that show the morphol- 
ogy of the high- velocity blue lobe gas. From this figure. 



it is clear that at higher outflow velocities the emission 
extends farther from the source, compared to lower ve- 
locities. In addition, in panels b through d the emission 
shows an elongated ring or shell-like structure. This im- 
plies that the protostellar wind not only entrains ambient 
gas along the parabolic walls visible at most blue shifted 
velocities, but also at angles closer to the outflow axis. 
This general behavior, in both morphology and velocity 
distribution, is consistent with the expected properties 
of a molecular outflow entrained by a w ide-angle wind 
dShu et al.||1991[ |Arce fc Goodmanll2QQ2l ) . 

In order to further investigate the underlying proto- 
stellar wide-angle wind and entrainment mechanism in 
the blue lobe, we constructed a position- velocity {p — 
digram along the axis of the outflow, shown in Figure 
This figure shows an inclined parabolic structure, as ex- 
pected for an outflow entrained by a wide-angle wind 
with a non-zero inclinatio n of the outflow axis w ith re- 
spect to the plane of sky ( Lee et al.||2000| |200lD . The 
relatively high angular resolution and sensitivity of our 
ALMA data results in one of the most (if not the most) 
clear case where such a parabolic structure is detected 
in the p — v diagram of a molecular outflow lobe. The 
structure in the blue lobe extends to outflow velocities of 
up to ~ —30 km s~^. 

We follow the simple analytical model by |Lee et ah] 
( |2000[ ) to describe the blue lobe of the HH 46/47 molec- 
ula r outflow. In this model, which is based on the study 
by Li & Shu (1996), the molecular outflow is made up 



of the gas that is swept-up by a wide-angle wind that 
propagates into a flattened (toroid-like) dense core. The 
resulting molecular outflow can then be modeled by a ra- 
dially expanding parabolic shell with a velocity structure 
in which velocity increases with distance (i.e., a Hubble 
law velocity structure). The mod e l is pi ctorially repre- 
sented in Figure 21 of Lee et al. (2000), and in cylin- 
drical coordinates the morphology of tne shell can be 
represented by the equation z = CR^, while the veloc- 
ity components of the shell along the z and R axes are 
represented by v^ = VqZ^ vr — VqR. In these equations 
C and Vq are free parameters with units of arcsec"^ and 
km s~^ arcsec"^, which are constrained by the shape of 
the molecular outflow shell and the shape of the p — v 
diagram, respectively. 

A map of the high- velocity blue lobe parabolic struc- 
ture was produced by integrating the emission for outflow 
velocities between —26 and —4 km s~^, in order to avoid 
"contaminating" emission from the cloud and the binary 
component outflow (see Figure [3|. We rotated the map 
by —30° and flt a parabola {z = Ci?^) to the structure. 
The flt gives a value for C of 0.3 ± 0.05, where the errors 
indicate the range of values that provide a reasonable flt 
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Fig. 4. — Velocity range-integrated intensity maps of the CO blueshifted emission. The velocity interval of integration is given on the 
top of each panel. In all panels the lowest contour and contour steps are 0.07 Jy beam"-*^ km s"-*^. The synthesized beam is shown on the 
lower left corner of each panel. The red cross shows the position of HH 47 IRS. 



to the map. We then fit the p — y diagr am with the wide- 
angle wind model of 'Lee et alj ( 2000 ) described above. 
We constrain the inclination angle with respect to the 
plane of the sky (i) to 30 ± 10°, close to the values de- 
rived by Eisloeffel fc Mundt^ (1994) and Hartigan et al.^ 

(2005), assuming a distance to the parent globule of 450 

1 



pc. We find that a value of Vq = 2.3 ± 0.2 km s " arc 
sec~^ and an value of z = 29 ± 1° result in a reasonable 
fit to the p — V cut. The errors indicate the range of 
values for which reasonable fits are attained, and choos- 
ing significantly different values for these two parameters 
would result in a parabolic curve with a different width 
or inclination compared to the p — v diagram obtained 
from our data. 

This simple analytical model is probably not the only 
kinematical model that may reproduce the observed ve- 
locity distribution. However, the fact that we can rea- 
sonably fit the integrated intensity map and the p — v di- 
agram with the model described above suggest that the 
observed blue lobe could very likely have been formed by 
the entrainment of a wide-angle wind and the ambient 
cloud. Our data, however, cannot distinguish between 
possible wind-launching mechanisms that can produce 
wide-angle winds (e.g., X-wind, Shang et al. 2006 or disk 
winds, Pudritz et al. 2006). We would need higher an- 
gular resolution data to probe down to less than ~ 50 
AU from the source in order to possibly discr iminate be- 
tween different launching mechanisms (e.g. 
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Fig. 5. — Position- velocity diagram along the axis of the blue 
lobe. The figure was constructed using a map with velocity resolu- 
tion of 0.5 km s"-*^, and summing the emission over a 7.5'' (about 
3 beams) wide cut, with a position angle of 60°. Contours have 
value of 0.35, 0.66, 1, 1.5, 2, 3 Jy beam"-*^, and subsequent steps of 
1 until 10 Jy beam~^, and thereon in steps of 2 Jy beam "-*^. The 
light blue parabola shows the model fit to the data (see §|4.1|). 



al. 2006). 



t might seem conflicting that the blue outflow lobe ap- 
pears to be entrained by a wide-angle wind, when it is co- 
incident with the base of a well-studied jet (see Figure [6|. 
However, it is possible for both a collimated wind compo- 
nent (i.e., a jet) to co-exist with a wide-angle wind 



g., m odels by [Shang et aE] 2Q07| Pudritz ( 
endt||2009[[Tomida et al ||2010p . in fact, thi 
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et aLllMyf 
is is not the 
tirst source to exhibit such characteristics. Three exam- 
ples where a "dual-component" wind have been invoked 
to explain the mm and optical/IR outflow observations 
are HH 111 (Nagar et al. 1997, Lee et al. 2000, however 
see Leflo ch et al. 2007 for an al ternative explan ation), 
HH 315 dA^e fc Good man||2002 '), and B5-IRS 1 ( [YtTet 
[aL]|1999 ). In the first two cases the protostellar source 
is close to the cloud edge, the blue optical jet (or HH 
flow) is clearly seen to reside outside the dense parts of 
the cloud, and the morphology of the molecular outflow's 
blue lobe is consistent with it being formed by a wide- 
angle wind (similar to HH 46/47). This is possible if 
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Fig. 6. — Comparison of blue CO lobe and optical image of HH 
46/47 jet. The HST i mage was tak en with the [SH] (F673N) fil- 
ter in 2008 by Hartig an et al.| ( |2011| . Contours show the emission 
integrated over the velocity range from -6.8 to -3.2 km s"-*^. This 
velocity range was chosen to avoid contamination from cloud emis- 
sion and BO. The lowest contour and subsequent contour steps are 
0.15 and 0.2 Jy beam"-*^ km s"-*^, respectively. The yellow cross 
shows the position of HH 47 IRS. The two dotted lines connect the 
position of the source to outermost positions of the wiggling jet, 
with respect to the outfiow axis. The dashed black line perpendic- 
ular to th e ou tfiow cavity shows the position of the p — v cut shown 
in Figure [Totx. 
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Fig. 7. — Velocity range-integrated intensity maps of the CO redshfited emission. The velocity interval of integration is given at the 
top of each panel. In the left, middle and right panels the lowest contour and contour steps are 0.2, 0.105, and 0.2 Jy beam"-*^ km s"-*^, 
respectively. The synthesized beam is shown on the lower left corner of each panel. The red cross shows the position oif HH 47 IRS. 



the underlying protostellar wind has both a cohimated 
(jet) and a wide-angle component and the gas in the 
blue lobe is mostly entrained by the wide-angle compo- 
nent. The reason that the wide-angle wind dominates the 
gas entrainment in these sources (including HH 46/47) 
is possibly due to the fact that the jets lie in a region of 
very low density molecular gas, on the outskirts of the 
cloud. Numerical simulations of jet (shock-driven) out- 
flows show that most of the gas entrainment takes place 
at (or near) the head of the bow-shock (e.g., Smith et| 



al 



1997 



Lee et al. 



2QQ1| ), while in a radially expand- 
ing wide-angle wind entrainment mostly takes place in a 
wide-angle shell originating at the source. Hence, in a 
dual component wind where the jet bow shocks mostly 
reside outside the cloud, the circumstellar molecular gas 
will mostly be accelerated by the wide-angle wind compo- 
nent close to the source, where there is enough molecular 
material for the wind to entrain and form the observed 
molecular outflow. 

Further evidence that the blue lobe of the HH 46/47 
molecular outflow is entrained by a dual-component 
wind is observed by comparing the morphology of the 
blueshifted gas and that of the optical jet. In Figure [6] we 
show contours of the low- velocity blueshifted lobe plotted 
over the HST image of the HH 47/46 jet. The northern 
wall of the CO blue lobe follows the northern edge of the 
optical nebula, which implies that the bright parabolic 
structure observed in CO traces the walls of the outflow 
cavity. The jet's wiggling structure suggests that t he jet 
axis changes with time, or precesses (i.e, Reipurth et al. 
|2000). In Figure [6] we show the opening angle of the pre- 
cession cone in trie plane of the sky, delimited by lines 
that connect HH 47 IRS and the bright emission knots at 
the most extreme angles with respect to the protostellar 
source. It is clear that the opening angle of the molec- 
ular outflow lobe is substantially wider than that of the 
jet precession cone. It is therefore highly unlikely that 
precession of the jet alone could produce the wide-angle 
cavity traced by the CO outflow. 

4.2. Red lobe 

Compared to the blue lobe, the red lobe of the HH 
46/47 molecular outflow exhibits a more complex spatial 
and kinematical structure. In Figure [7| we show maps of 
the redshifted outflow emission integrated over different 
velocity ranges chosen to highlight the different impor- 



tant structures discernible at different velocities. These 
are discussed in detail below. 

4.2.1. Morphological evidence for the existence of a 
wide-angle wind 

Figure |7|a presents a map of the low outflow velocity 
emission (1.2 < Vout < 2.7 km s~^), where the lobe shows 
a clear wide-angle structure, with an opening angle of 
65°, within 40" of the source. This wide-angle structure 
is drastically different from the morphology of the bright 
IR ou tflow cavity walls seen in the IRAC Spitzer observa- 
tions ( Noriega-Crespo et al.|2QQ"4 ), which is similar to the 
structure of the CO outflow red lobe at distances from 
the source further than 40" (see Figure Isl). The loop-like 
structure seen in the IRAC images has oeen adequately 
modeled as the c avity walls of a jet-driven outflow by 
Raga et al. (2004). This model, however, does not repro- 
duce the wide-angle structure seen in low- velocity CO 
emission nor the wide-angle IR nebula seen close to the 
source in re-processed IRAC images of HH 46/47 (see 
below). It seems, therefore, that the wide-angle struc- 
ture is not produced by the collimated (jet) component 
responsible for the bright outflow cavity walls seen both 




Right Ascension (J2000) 

Fig. 8. — Comparison of red CO lobe with Spitzer IRAC 2 
(4.5 /im) image. The Spitzer data is from iNoriega-Crespo et al.] 
(1^04) and has been reprocessed with a deconvbrut ion algorithm 
to reach an angular resolution of ~ 0.6 — 0.8^^, with 60 iterations 
(see Noriega-Crespo &; Raga 2012 for details). White contours 
show the integrated intensity of the CO redshifted lobe (the same 
as shown in Figure IT]). Black contours show the 3 mm continuum 
emission around Hir47 IRS. The lowest contour and subsequent 
contour steps are 2 mJy beam"-'^. The white dash line shows the 
position of the p — v cut shown in Figure [Tot ). 
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in CO and IR, extending form about 30'^ to 2' away from 
the source. Similar to the blue lobe, the red lobe shows 
evidence of both a collimated and a wide-angle wind com- 
ponent. 

Recently, Velusamy et al. (2007) presented reprocessed 
images of tne IK AC observations of HH 46/47, that in- 
cluded suppression of side lobes from bright sources and 
enhanced resolution (compared to the images presented 
by Noriega- Crespo et al. 2004). The 3.6 and 4.5 /im 
images show diffuse and extended emission close to the 
source with a parabolic morphology, with an op ening an- 
gle of about 110°, which Velusamy et al. (2007|) argue is 
due to scattered light arising from a wide-angle outflow 
cavity (see Figure |8|. This structure is wider than the 
CO outflow lobe, and we suspect we are not able to trace 
the full extent of the IR reflection nebula due to the high 
opacity of the CO (1-0) at velocities close the ambient 
cloud. In the few cases where a molecular outflow has 
been mapped in ^^CO, it typically shows a wider opening 
angle (at lower outflow velocities) than the lobe traced 
by the ^^CO emission and the optical/IR reflection neb 
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ula (i.e,IO hashi et al.|1997{ 'Tafalia 

|Sargent[5006) . We expect that future observations using 
lines with a lower opacity (or that probe higher densi- 
ties) will be able to trace the wide angle cavity outlined 
by the faint extended IR emission. 

4.2.2. Evidence for episodic mass ejections 

Among the most distinct structures in the velocity- 
integrated map shown in Figure^ are three bright and 
compact features along the outflow axis. These local- 
maxima of emission are the same clumps identified in 
Figure [l] as Rl, R2, and R3. They are more clearly seen 
in Figure [7^ because of the limited velocity range used to 
produce this figure compared to Figure [l] Table [^ shows 
the distance from the source as well as other physical 
properties of each of the molecular clumps. The distance 
from the source to Rl (about 0.1 pc) is approximately 
equal to the distance between Rl and R2, and the dis- 
tance between R2 and R3. This suggests they might be 
produced by periodic ejection episodes, similar to those 
found in other molecular outfiows (e.g., L1448, Bachiller 
et al. 90; RNO 43 Bence et al. 96; HH 300, Arce & Good- 
man 2001b). In fact, Rl and R2 are at similar distances 
from HH 47 IRS as the bright bow shock structures HH 
47B and HH 47A, respectively, along the optical jet, in 
the opposite (blue) lobe of HH 46/47. These two HH 
knots are thought to arise from differ ent mass ejection 
episodes ( [Reipurth fc Heathcote 1991| ), and the Rl and 
R2 may be associated with the counter-ejections (in the 
redshifted lobe) of the mass-loss events that caused HH 
47B and HH 47A. Similarly, R3 is coincident with HH 
47C (in the red lobe), which is thought to be associated 
with the ejection responsible for the HH 47D bow shock 
(in the blue lobe). 

To further investigate the nature of the outfiow clumps 
we show, in Figure [9) a p — v diagram along the axis 
of the redshifted molecular outfiow lobe. At distances 
greater than 30^' from the source, the p — v diagram 
shows three regions where there is a clear increase in the 
outfiow velocity coincident with the position of the red 
outfiow clumps. These features in the p — v diagram, gen- 
erally referred to as "Hubble wedges" ( Arce & Goodman 
|2001a) , show maximum outfiow velocity increasing as a 



TABLE 2 

Properties of clumps in redshifted CO 

LOBE 


OUTFLOW 


Clump 


Distance from source^ 

(arcsec) (pc) 


Age^ 

(yrs) 


Mass^^ 
(10-3 Mo) 


Rl 
R2 
R3 


42 0.1 
81 0.2 
118 0.3 


360 
650 
945 


0.9 
0.8 
0.3 



^ Values shown have been corrected for the inclination of the 
outflow axis, assuming i = 30° . 

Mass estimate obtained assuming optically thin emission 
and an excitation temperature of 15 K. Values shown should 
be treated as lower limits. 

function of distance from the source, within a localized 
length along the outfiow axis. The general consensus is 
that such features are formed by a collimated (jet-like) 
episodic fiow with a significantly varying mass-ejection 
rate, where each Hubble wedge is caused by the bow- 
shock entrainment of an individual mass -ejection event 
( |Arce fc Goodman|[2QQla||Lee et aL]|2QQlD . We therefore 
interpret the redshifted outflow clumps m HH 46/47 as 
arising from different mass ejection episodes from HH 47 
IRS. 

In Table [2] we list estimates of the timescale (i.e., age) 
for each mass ejection episode in HH 46/47 assuming that 
all clumps are moving away from the source at a constant 
space velocity. We expect the ejected material will de- 
celerate as it moves through the cloud and interacts with 
the surrounding medium, which will res ult in in slower 
veloc i ties for older ejected ma terial (e.g., Cabrit & Raga 
2QQQ[ Goodman fc Arce||2QQ4| ). Yet, we can use an aver- 
age now velocity of the jet [Vjet) for all outfiow clumps to 
obtain an estimate o the ejection age for each mass ejec- 
tion episode. Existing studies indicate that for HH 46/47 
Viet ^ 30 km s~^ (Eisloe ffel fc Mundt ^ ^l994( [Micono et 
al. 1998] [Hartigan et al. ' 7UU5[ ). Usmg this velocity, we 
estimate that a major miass ejection event in HH 46/47 
takes place approximately every 300 yrs, consistent with 
the results from studies of the optical (blueshifted) jet. 

The CO outfiow clumps we detect are most likely made 
of gas that has been entrained near the current posi- 
tion of the clump. Recent observations indicate that the 
outfiow mass r ate in HH 46/47 is, at the most, abou t 



10"_MrD yr" (Antoniucci et al. 



2008 



Garcia Lopez et 



aL||2010 |). From the length of K2 along the outflow axis 
(i.e., ^ 20'^), and assuming a constant Vj^t = 300 km s""*^, 
we estimate that the episodes of high mass ejection are 
about 10^ yrs. With the mass outfiow rate given above 
this would result in a total of 10~^ Mq ejected in the 
wind, which is about two orders of magnitude less than 
the estimated gas mass of R2. The mass loss rate during 
the episode associated with R2 would have to be about 
two orders of magnitude higher than the outfiow mass 
loss rate in the most recent mass ejection episode traced 
by the IR jet within 10^' of the source (i.e., the IR knots 
Zl and Z2 in the red lobe, where Garcia Lopez et al. 2010 
measure the mass outfiow rate), for the wind mass to be 
similar to the estimated gas mass of R2. Hence, the vast 
majority of the detected molecular gas in the redshifted 
outfiow clumps is cloud material that has been entrained 
by the protostellar wind (and not material from the wind 
itself). 
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Fig. 9. — Position- velocity diagram along the axis of the red lobe. 
The figure was constructed using a map with velocity resolution 
of 0.08 km s"-*^, with a position angle of 60°, and summing the 
emission over a 7.5^^-wide cut. The lowest contour and subsequent 
contour steps are 1 and 1.5 Jy beam"-*^, respectively. The positions 
of the redshifted outflow clumps along the outflow axis are shown. 



Episodes of enhanced mass outflow rate are thought 
to be caused by an increase in the mass accretion rate 
onto the protostar. In the case of HH 46/47, there seems 
to be a mechanism that is able to produce an enhance- 
ment in the accretion rate every few hundred years. This 
timescale is significantly shorter than the estimated or- 
bital period of the binaries in the HH 47 IRS system 
(see Reipurth et al.p QQQ). This rules out triggering of 
gravitational instabilities in the circumstellar disk of the 
source that drives the HH 46/47 outfiow by the close 
passage of the known binary companion as a mechanism 
for enhanced disk accretion. A closer (yet unobserved) 
companio n would have t o exist for companion-disk inter- 
actions (Reipurth"2000') to be a feasible mechanism for 
triggering episodic mass ejections. Other mechanisms 
that could produce episodic accretion (and hence mass 
ejection) events have been discussed in the literature, 
where disk instabiliti es are the main cause of enhanced 
accretion events (e.g.,'Vor Qbyov fc Basu||2QQ5 Zhu et al._ 
|2010; Machida et al. 201ll). In many ot the models the' 
timescales between successive episode may be very dif- 
ferent depending on the assumed input parameters, and 
in some of these models successive events of high mass 
accretion rate can occur on timescales as low as a few 100 
yrs to about 10^ yr, consistent with the observed proper- 
ties of the HH 46/47 outfiow. Future studies of outfiow 
episodicity using a large sample may help constrain vari- 
able accretion models. 

Another morphological feature detected in the red- 
shifted molecular outfiow lobe that may be caused by 
multiple mass ejections is the bifurcation of the northern 
cavity wall. This feature is seen at low and medium out- 
fiow velocities (Figures [Tja and[7|)) at about V away from 
the protostar along the northern wall of the outfiow lobe. 
Here the wall appears to divide into two filaments: a 
bright component that extends southwest approximately 
parallel to the outfiow axis (PA ~ —120°) and a fainter 
component that makes a sharp southward turn towards 
the position of R2 (see FigureJTJa and b). Models and sim- 
ulations of molecular outfiows formed by episodic winds 
show that a thin shell or cavity is formed around each 
i nternal bow shock formed by differe nt ejection episodes 
(jGueth et al.| |1996| |Lee et al.| |2QQT] ). We thus surmise 
that the filamentary feature that connects R2 and the 
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Fig. 10. — Position-velocity diagrams perpendicular to the out- 
flow axis for the blue (left panel) and red (right panel) lobes. In 
both p — V diagrams the offsets are given with respect to the posi- 
tion of the outflow axis in the cut. The location of the p — v cuts for 
the blue and red lobes are shown in Figures [6] and [s] respectively. 
In the left panel the lowest contour and the subsequent contour 
steps are 0.3 and 0.15 Jy beam"-*^, respectively. In the right panel 
the lowest contour and the subsequent contour steps are 0.07 and 
0.16 Jy beam"-*^, respectively. 

northern wall of the red lobe traces the shell formed by 
the interaction of the ejection episode associated with R2 
and the ambient gas. Faint emission connecting Rl to 
the southern wall of the lobe may indicate that a similar 
structure is also driven by Rl. 

4.2.3. High-velocity redshifted emission 

Figure [7|c shows a velocity-integrated map over high 
CO outfiow velocities in the red lobe (6.2 < Vout < 
21.2 km s~^). Here the outfiow opening angle near the 
source (~ 50°) and the lobe width are smaller than at low 
outfiow velocities. The overall arc-like morphology of the 
northern wall is consistent with the shape expected for 
an o utfiow lobe formed by (jet) bow-shoc k entrainment 
(e.g.,|Raga fc Cabrit | 1993| | Lee et al. | 200lD . This and the 
other morphological and kinematical features discussed 
above imply that the dominant entrainment mechanism 
of the CO outfiow observed beyond 40^' arises from a col- 
limated episodic wind with a series of bow shocks along 
the axis. On the other hand, near the source, and at low 
outfiow velocities the molecular outfiow appears to have 
been formed by the entrainment from a wide-angle wind. 
A possible explanation is that the wide-angle component 
has a lower velocity (and possibly a lower density) than 
the collimated (jet-like and episodic) component. Near 
the source, where the ambient density is high, both the 
wide-angle and collimated component will entrain the 
gas, but the effects of the wide-angle wind will be seen 
most clearly at low outfiow velocities and at relatively 
large angles away from the outfiow axis. The faster mov- 
ing (and denser) collimated component is able to travel 
further away from the source (and even puncture out of 
cloud near the position of R3), and thus the CO outfiow 
further out from the source (and at high-velocities) will 
exhibit properties of it being mostly entrained by a col- 
limated wind. Certainly, a successful model for the HH 
46/47 molecular outfiow would need to reproduce the dif- 
ference in morphology between the low and high outfiow 
velocities and at different distances from the source. 
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4.3. Velocity Gradients in Outflow Cavity 

The high sensitivity, high velocity resolution data al- 
lowed us to study the kinematics of the outflow cavity in 
detail. In Figure [lO] we plot the p — v diagram made by 
cuts perpendicular to the outflow cavities of the blue and 
red lobes. In the p — v diagram for both cavities there is 
a wide range of velocities present at the position of the 
cavity walls, and a gradual velocity gradient with most 
emission near the outflow axis at higher velocities than 
(most of) the emission closer to the outflow walls. In 
the blue lobe, at the position of the outflow cavity wall, 
there is detectable CO at velocities from Vout ^ — 5 km 
s~^ to Vout ^ —0.5 km s~^. Yet, most of the emission at 
this position is at velocities close to that of the ambient 
cloud (i.e. , between Vout of —1.5 and —0.5 km s~^, see 



regions near the cavity walls (e.g., Arce fc Sargent|2QQ4). 



Figure 10 a). This in contrast to positions closer to the 



outflow axis, where the bulk of the emission is at slightly 
higher velocities, between Vout of —2.0 and —3.0 km s~^. 
A similar velocity distribution, but with much faster ve- 
locities, has been observed in the HH 46/47 optical jet 
with fast radial velocities (~ —100 to —150 km s~^) in 
the jet core and s lower velocities (^ —25 to —50 km 
s"^) at the edges ( [Morse et aLl|l994| ). In the CO out- 
flow red lobe we see a spread in velocities (i.e., emission 
over a range of Vout of about 2 km s""*^) at the position 
of the outer (wide-angle) cavity and the inner cavity, as 
well as an increase in the average (and minimum) ve- 
locity of the CO with distance from the outer outflow 
cavity wall: from velocities cloud to that of the ambi- 
ent cloud {vout^ 0-5 km s~^) at the outer cavity wall, 
to Vput^ 1.6 km s~^ at center of the outflow lobe (see 



Figure [10^). 

The wide spread in velocities at the position of the cav- 
ity wall may be explained by both geometric effects and 
the expected interaction between the outflow and the 
surrounding cloud. Observations with limited angular 
resolution are expected to detect a large range in radial 
velocities (i.e., the component of the gas velocity along 
the line-of-sight) at the limb-brightened edge of an ex- 
panding parabolic (or cylindrical) cavity, as outflow gas 
from the front and back side of the cavity wall should 
contribute to the total detected emission. In addition, 
given that the cavity walls is where the outflow-cloud 
interaction takes place, one would expect the gas to be 
more turbulent than close to the cavity center, and hence 
show a wider spread in velocities. 

It is not clear what causes the observed velocity gra- 
dient perpendicular to the outflow axis. The decrease in 
velocity of about 1 km s~^ from the center of the cav- 
ity to the cavity wall may be due to different processes. 
One possibility is that this behavior in the CO outflow 
results from the entrainment of a protostellar wind that 
has higher velocities towards the wind axis. For example, 
a disk wind where the velocity of the wind depend s on 
the radius at which is launched ( Pudritz et al.|200'6 ) or a 
jet made of many small bow-shocks where the entrained 
material at the head of the bow shocks (along the center 
of the jet) is expected to have higher velocities than the 
material entrained at the side wings of the bow shock 
( [Heathcote et al.||1996 ). Another alternative is that the 
velocity gradient results from the momentum-conserving 
wind entrainment that produces the molecular outflow, 
where lower CO outflow velocities are expected in denser 



An additional possibility is that this is an ettect due to 
the geometry of an expanding parabolic, conical or cylin- 
drical cavity inclined with respect to the plane of the 
sky, where CO emission with larger components of the 
velocity along the line of sight are preferentially detected 
closer to the center of the cavity. Undoubtedly, the de- 
tected gradient is not due to outflow rotation, as it has 
be en recently claimed for a few other sources (e.g., Lee et 



■al.l2009J [Zapata et ar]|2QTQl [Choi et al.[[2QTT| |Fech et al. 
2012). The observed decrease in velocity from the out- 
tiow axis towards the cavity walls is in contrast with the 
expected observational signatures of outflow rotation: a 
gradient across the entire lobe length with one side of the 
lobe showing bluer (or redder) velocities with respect to 
the other side of the lobe (e.g., [Launhardt et al.[[2009 ). 



4.4. Outflow from binary 

As described in § [3.2. 1[ the channel maps show a pro- 
tuberance extending about 20^' southeast of the source, 
at low blueshifted outflow velocities {vout ^ —1.8 km 
s~^). We believe that this clump (which we name BO) 
is not related to the main blueshifted outflow lobe, as it 
does not follow the general parabolic shape seen at most 
blue shifted velocities. A p — v diagram oriented along a 
line from HH 47 IRS to BO, at PA - 132° (Figure [n]), 
shows that the maximum velocity increases from near 
the source to the position of BO, similar to the Hubble- 
wedges of the red outflow clumps. We speculate that BO 
traces a region where an underlying protostellar wind is 
entraining molecular gas. 

A recent study of the HH 46/47 jet shows that the 
proper motion on the plane of the sky of one of the most 
recent ejections (within ~ 3'^ of HH 47 IRS) is along a 
line with a position angle of about 40° on the plane of 
the sky, while further away (~ 35^0 from the source th e 
jet motion has a PA closer to 55° (Harti gan et al.[[2011 ). 
Thus, it is very unlikely that the ejection responsible lor 
BO is part of the flow responsible for the HH 46/47 jet, 
and the most likely candidate for powering the outflow 
that drives BO is the other binary companion in the HH 
47 IRS binary system not responsible for the main HH 
46/47 outflow. Weaker and smaller outflows from binary 
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Fig. 11. — Position- velocity diagram along the axis of the pro- 
posed outflow powered by the binary companion in the HH 47 IRS 
system not responsible for the main HH 46/47 outflow. The lo- 
cation of the p — v cut is shown in Figure ^ The lowest contour 
and the subsequent contour steps are 0.1 and 0.3 Jy beam"-*^, re- 
spectively. The locations of the blueshifted outflow clump BO and 
clump A, along the p — v cut, are shown. 
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components have been observed in other sources (e.g., 
HH 111, Cernicharo & Reipurth 1996; BHR 71, Parise 
et al. 2006; IRAM 04191+1522, Lee et al. 2005, Chen et 
al. 2012), and we believe that HH 46/47 shows another 
case of drastically different outflows arising from bina- 
ries in the same system. Current observations of this 
system do not have the required angular resolution to 
resolve which outflow is powered by each of the binary 
components in HH 47 IRS. 

Figure 11 shows that there is a clump (which we name 
A) with very low red (and blue) outflow velocities at 
about 30'^ northwest of HH 47 IRS on what would be 
the counter lobe of the outflow associated with BO. The 
distance from the source to A is slightly more than the 
distance to BO, and the position of A is coincident with 
the edge of the globule seen in our CO channel maps 
(Figure Isl). It is not clear if clump A arises from the 
counter ejection of BO or if it is just cloud emission. Fur- 
ther observations are needed to ascertain the nature of 
this clump. 

5. SUMMARY & CONCLUSIONS 

We present the first interferometric map of the molec- 
ular (CO) outflow associated with the well-known HH 
46/47 flow, powered by HH 47 IRS. With our high veloc- 
ity (0.08 km s"^) and angular (~ 3'^) resolution ALMA 
Cycle observations we have investigated the kinematics 
and morphology of the molecular outflow. The 4.2' long 
mosaic map of the CO (1-0) emission covers the length 
of the HH 46/47 outflow from HH 47 A, in the northeast, 
to HH 47C, in the southwest, and provides the data to 
study both the blue- and redshifted molecular outflow 
lobes in detail. Our main conclusions are as follows. 

The unprecedented sensitivity of our data enable us to 
detect outflow emission at much higher velocities than 
previous CO(l-O) studies of the source; up to outflow ve- 
locities of about of —30 and 40 km s~^ in the blue and red 
lobes, respectively. Even though the very high velocity 
(I'^on^l > 20 km s~^) emission is constrained to about 15'^ 
(or 7000 AU) from the source, detection of this outflow 
emission results in significantly higher values of the out- 
flow kinetic energy and momentum, compared to results 
from existing single dish CO (1-0) observations. If HH 
46/47 is representative of molecular outflows from low- 
mass stars, our results imply that other similar molecu- 
lar outflows may be much more energetic and carry more 
momentum than previously thought. The HH 46/47 flow 
is the first molecular outflow from a low-mass star to be 
mapped in the ALMA Cycle phase and subsequent 
ALMA observations using more antennas will provide 
even more sensitive maps. From our results we expect 
that future outflow studies with ALMA will show that 
outflows can have much more impact on their surround- 
ing cloud than previously thought. 

The blue and red molecular outflow lobes show very 
distinct morphologies and kinematics. The blueshifted 
lobe only extends up to about 30'' (14000 AU) north- 
east from HH 47 IRS. The compact size of the molecular 
outflow lobe is due to the fact that the outflow-powering 
protostar resides at the edge of its parent globule and the 
protostellar wind is only able to push and accelerate a 
limited amount of molecular material before it leaves the 
cloud. The blueshifted lobe shows a very clear parabolic 
structure in both the integrated intensity map and in the 



position- velocity diagram along the outflow axis. The 
morphology and velocity distribution are consistent with 
models where the molecular outflow is formed by the en- 
trainment of an underlying wide-angle protostellar wind. 
We argue that the co-existence of a wide-angle molecular 
outflow and an optical jet in the blue lobe of HH 46/47 
is possible if the underlying protostellar wind has both 
a collimated (jet) and a wide-angle component. In the 
blue lobe, the wide-angle wind component dominates the 
gas entrainment because the jet bow shocks lie in a re- 
gion with very low ambient density while the wide-angle 
component interacts with the the environment very close 
to the source, where there is enough ambient (molecular) 
gas for the wide-angle wind to entrain and form the ob- 
served molecular outflow. 

The redshifted lobe is about 2' (~ 0.3 pc) in length, ex- 
tending out to the southwest edge of the globule. Within 
approximately 40" (18000 AU) of the source the red lobe 
shows a wide-angle structure which is very different from 
the loop-like structure seen in IR images and the CO lobe 
structure at greater distances from the source. Our re- 
sults also show that, beyond 40" and along the axis of the 
red lobe, the HH 46/47 outflow has three clumps clearly 
detected in the integrated intensity map and position- 
velocity diagram. Their spatial distribution and velocity 
structure indicate that these clumps arise from prompt 
entrainment, most probably produced by bow shocks, 
arising from quasi-periodic episodes of increased mass 
ejection in HH 47 IRS. We argue that, similar to the 
blue lobe, the underlying protostellar wind has both a 
collimated and a wide-angle component and that close 
to the protostar the molecular outflow red lobe is dom- 
inantly driven by the wide-angle component while the 
CO outflow observed beyond 40" mainly arises from the 
entrainment by a collimated episodic wind. The wide- 
angle wind of this, or any other, outflow may affect the 
infall and the protostar's mass-assembling process, and 
future ALMA observations with higher angular resolu- 
tion and using higher density tracers should concentrate 
on studying the impact of the wide-angle wind on the 
inf ailing core and circumstellar envelope. 

The position of the three molecular outflow clumps 
along the axis of the redshifted lobe suggests that they 
are caused by the same ejection event that produced 
three major bow shocks in the blue lobe observed in 
the optical. Our results confirm previous studies that 
indicate that HH 47 IRS goes through episodes of in- 
creased mass outflow (and accretion) rate every few hun- 
dred years. Future (higher resolution) multi-epoch ob- 
servations with ALMA may be able to detect the proper 
motion of the bow shock-driven outflow clumps in HH 
46/47 and other molecular outflows. These studies could 
help in our understanding of outflow variability and con- 
strain variable accretion models. 

Our high velocity resolution data allowed us to study 
the kinematics of the outflow cavity and search for sig- 
natures of outflow rotation. We find, in both blue and 
red lobes, there is outflow emission over a wide range 
of velocities at the cavity walls, and the existence of a 
gradual gradient where the velocity increases towards the 
outflow axis. It is not clear what produces this velocity 
structure. Yet, we are certain it is not consistent with 
outflow rotation. Future comparison of our data with 
numerical models will allow us to investigate the origins 
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of this kinematic behavior in the entrained gas. 

Our CO (1-0) data detect a smah clump at blueshifted 
velocities southeast of HH 47 IRS, which we believe arises 
from the interaction of a protostellar wind driven by the 
binary component in HH 47 IRS that does not drive the 
main HH 46/47 outflow. The clump, which lies at a 
position angle very different from that of the main HH 
46/47 flow, has a velocity distribution consistent with it 
being formed by bow shock prompt entrainment. Future 
observations with an order of magnitude higher angular 
resolution than the maps presented here will be able to 
confirm our results and resolve which outflow is powered 
by each of the binary components in HH 47 IRS. 
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